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The study of 7-ray emission from galactic sources such as supernova remnants (SNR) may provide key insights 
into their potential role as accelerators of cosmic rays up to the knee ( ~ 10^^ eV). The VERITAS Observatory 
is sensitive to galactic and extragalactic 7-ray sources in the 100 GeV to 30 TeV energy range. We report here 
on VERITAS observations of the vicinity of the cocoon of freshly accelerated cosmic rays reported by Fermi, 
which lies between potential accelerators in the Cygnus OB2 association and the 7-Cygni SNR. A particular 
focus is placed on the source VER J2019 +407 in 7-Cygni. 



1. The VERITAS Instrument 

VERITAS, located at the Fred Lawrence Whipple 
Observatory near Tucson, Arizona, is an array of four 
12-meter imaging atmospheric Cherenkov telescopes. 
Each telescope has a pixelated camera comprised of 
499 photomultiplier tubes with a 3.5° field of view. 
Designed to detect photons of astrophysical origin be- 
tween 100 GeV and 30 TeV, VERITAS detects and 
images the secondary Cherenkov light produced when 
gamma rays and cosmic rays initiate particle cascades 
in the upper atmosphere. Stereoscopic reconstruction 
of events using multiple telescopes allows for single- 
photon angular resolution of better than 0.1° and en- 
ergy resolution on the order of 15-25%. In its cur- 
rent configuration, VERITAS can detect a source with 
~ 1% of the Crab Nebula fiux in less than 30 hours. 
The instrument has been operating in full array mode 
since 2007. See [Holder et aLll2006i for further details 
regarding the operation of VERITAS. 



2. The 7-Cygni Supernova Remnant 

SNR G78. 2-^2.1, also known as the 7-Cygni 
supernova remnant (SNR), is a ^ 1° diameter 
shell-like radio and X-ray SNR[Higgs et all Il977l 
iLozinskava et aLll2000ll. It is estima t ed to be at a dis- 
tance of ~ 1.7 kpc fHiggs et al.|[T977l . ILozinskava et al.l 
|200C^ an d to b e approximatel y 5000 — 7000 years 
old^Hi ggs et al.l [l97 7. Land ccker et al.l [1980] . It is 
thought to be in an early p hase of adiabatic expans ion 
into a low-density medium Lozinskava et al.ll2000j |. A 
slowly expanding H i shell im mediately surround - 
ing the radio shell was found bv lGosachinskiil (200 1| . 
which ILozinskava et al.1 j2000l | believes to have been 
created by the progenitor stellar wind. 

The radio she ll divides roughly into northern and 
south ern arcs Uchivama et al.l l2002l IZhang et al.l 
I1997 1 . Enhanced ther mal X-ray emission in the north 



suggests shocked gas [Uchivama et a l. 2002] and falls 
in a void of CO emission [Ladouceur fc Pi ncault 2008] . 



satellite Fermi, operating at GeV energies, has discov- 
ered a 7-ray pulsar PSR J2021-j -4026 at the center of 
the remnant f Abdo et al.ll2010al lb|. While this pulsar 
has a low luminosity (1.1 x 10^^ erg s~^) and a spin- 
down age (76.8 kyr) much greater than the estimated 
age of SNR G78.2-t-2.1, the pulsar kinematics make it 
probable that PSR J202 1-1-4026 was born with some- 
thing close to its current spin period and is the rem- 
nant of SNR G78.2-h2.1's progenitor star jTrepl et al.l 
l2010f . Diffuse 7-ray emission above 10 GeV is also 
repor ted by Fermi over the full extent of the rem- 



nant 



Lande et al. I l2012j . A point source co-located 



with VER J2019-I-407 was previo usly reported in the 
first and second Fermi catalo gs lAbdo et al.l l2010al 
iNolan et al.ll2012| ]. but , Lande et al.l j2012| c"onclude it 
to be artifact. 



3. The Cygnus Cocoon 



Strong optical emi ssion with sulfur lines also charac- 
terizes the region [MavromatakisI [2003| ] . The 7-ray 



lAckermann et all [2011 reported an extended re- 
gion of emission above a few GeV that they inter- 
preted as a freshly accelerated cocoon of cosmic rays. 
The cavity defining the cocoon region is outlined by 
ionization fronts visible in the mid- infrared. 7-Cygni 
overlaps one end of the cocoon and is a potential 
source of the cocoon's trapped cosmic rays. 



4. Detection of 7-ray Emission from the 
Direction of 7-Cygni 

Figure [1] displays the acceptance-corrected very- 
high-energy (VHE) 7-ra y excess map of the area 
around SNR G78.2-h2.1 jAliu el al.l [2013| . An ex- 
tended source is seen overlapping the northern edge 
of the remnant, with a detection significance of 7.5 
standard deviations. We use a binned, extended 
maximum-likelihood fit to the raw counts map to as- 
sess the VER J20 19-1-407 morphology. The source is 
modeled in the fit as a symmetric, two-dimensional 
Gaussian convolved with the VER J2019-I-407 point- 
spread function (PSF); the background is assumed 
to be flat before exposure effects are taken into 
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Figure 1: VERITAS 7-ray image of SNR G78. 2+2.1 showing the detection of VER J2019+407 and its fitted extent 
(black dashed circle). CGPS 1420 MHz co ntinuum radio contours at brightness temperatures of 23. 6K, 33. OK, 39. 6K, 
50K and lOOK (white) [Taylor et alJl20d3 ] outline the radio SNR; the star symbol shows the location of the central 
7-ray pulsar PSR J2021+4026. The inverted triangle and dot-dashed circle (yellow) show the fitted centroid and extent 
of the emission detected by Fermi above 10 GeV. The open and filled triangles (black) show the positions of Fermi 
catalog sources IFGL J2020.0-I-4049 and 2FGL J2019. 1-1-4040. The 0.16, 0.24, and 0.32 photons/bin contours of the 
Fermi detection of the Cygnus cocoon are shown in cyan. The white circle (bottom right corner) indicates the 68% 
containment size of the VERITAS 7-ray PSF for this analysis. 



account. Wc find a fitted extension of the two- 
dimensional Gaussian to be 0.23° ± 0.03°tjjt j;Q;g20sys, 
with fitted centroid coordinates R.A. 20''20™04.8^ 
Decl. -K40°45'36" (J2000). The statistical uncer- 
tainty in the centroid location is 0.03°, with a sys- 
tematic uncertainty of 0.018°. The systematic un- 
certainty considers both the telesco pe pointing erro r 
and systematic errors of the fit itself [Aliu el al.ll2013j | . 
The positions of the 7-ray pulsar PSR J202 1-^4026 
(IFGL J2021.5-h4026) (- 0.5° from VER J2019-H407 
) and the centroid of the emission above 10 GeV from 
the remnant (as seen by Fermi ) are also shown for 
reference. 

Figure O shows the spectrum of reconstructed 7- 
ray events within 0.24° from R.A. 20^19^48% Decl. 
-|-40°54'00". Runs where only three of four tele- 
scopes were operati onal h ave been excluded from this 
sample [Aliu el al.l l2013j | . The photon spectrum is 
consistent with a differential power law in energy. 



dN/dE = No X (E/TeV)"^, between the analysis 
threshold of 320 GeV and 10 TeV. The photon index 
is r = 2.37 ± 0.14stat ± 0.20sys and the flux normal- 
ization at 1 TeV is Nq = (1.5 ± 0.2stat ± 0.4sys) x 
10~^^ ph TeV~^ cm^^ s^^. The integral flux above 
320 GeV (5.2 ± 0.8stat ± 1.4sys x lO'^^ ^-^ ^^-2 g-i) 
corresponds to 3.7 % of the Crab Nebula flux above 
that energy (Aliu el al.ll2013} . 



4.1. 7-Cygni in X-Rays 

Figure [3] illustrates the region of en hanced X- 
ray emission overlapping VER J2019-t-407 [Aliu el all 
i2013i |. It displays the 0.7 — 3.0 keV exposure- 
corrected X-ray map using data from ASCA Se- 
quence #25010000 ( data originally presented by 
l Uchivama eTall j2002j |). generated by co-adding data 
from the two gas imaging spectrometers. A spectrum 
was extracted from a 12' x 24' elliptical region enclos- 
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Figure 2: Spectrum of VER J2019+407, derived from 
4-telescope data only. Points are the VERITAS 
spectrum; the line is a power-law fit with a spectral index 
of r = 2.37±0.14stat ±0.20sys and a flux normalization of 
No = (1.5 ± 0.2stat ± 0.4sys) X 10"^^ ph TeV"^ cm^^ s'K 
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Figure 3: ASCA X-ray view of G78. 2-1-2.1 between 1 
and 3 keV, overlaid with the VER J2019-f 407 smoothed 
photon excess contours (100, 150, 210 and 260 photons). 
The region used to extract a spectrum and the 
corresponding background region to the south of the 
remnant are indicated by white solid and dashed ellipses, 
respectively. A white star marks the position of 
PSR J202H-4026. 



5. Interpretation 



ing most of the X-ray emission inside the VERITAS 
contours, centered on coordinates R.A. 20^ 20™ IT'', 
Decl. +m° 45' 41" (J2000) and oriented with position 
angle 60°. We selected background photons from an 
identically sized ellipse near the center of the remnant 
at R.A. 20^^ 19™ 38% Decl. -^40° 27' 02" (J2000), with 
position angle is 130°. The source and background re- 
gions are displayed in Figure [31 

The source spec trum between .7 and 3.0 keV is 
shown in Figure l4l (Aliu el al.ll2013j . The X-ray spec- 
trum was modeled using an absorbed Raymond- Smith 
thermal plasma model with a best-fit temperature 
of kT = 0.57 ± 0.14 keV. With this model the col- 
umn density is Nh = (3.7 ± 2.0) x lO^^cm^^^ i^j^e 
normalization is N = 1.8 x 10^"^cm^^m and the 



-12 . 



absorption-corrected flux is 6.0 x 10^"' erg cm 's 
in the 0.5—8.0 keV band. This result differs sig - 
nificantly from that given in Uchivama et al.l [2002j | . 
which claimed an additional power-law component 
and a large Ne IX line feature. The divergence in re- 
sults appears to hinge, not on the choice of source re- 
gion, which is similar, b ut on the choice of background 
region |Aliu el al1l2013j | . We selected a background re- 
gion that was as close as possible to the source while 
Uchivama et al.|[2002j | chose a region 3.5° away. We 



can produce Uchiyama's results by choosing a back- 
ground r egion similar to U chiyama's. For further de- 
tails, see lAliu el all [2013 1. 



It is plausible that the VHE 7-ray emission seen 
from 7-Cygni arises from particles accelerated in 
shocks occuring at the interaction of the supernova 
ejecta and the surrounding medium. These parti- 
cles could be either accelerated electrons, which would 
produce the emission via inverse- Compton scattering, 
or accelerated nuclei. Should they be high-energy 
electrons, they would also be expected to produce 
X-ray synchrotron radiation, which would appear as 
a non-thermal power-law component in the X-ray 
spectrum. While our analysis of the ASCA X-ray 
spectrum does not argue for a non-thermal compo- 
nent, our upper limit on this component is still weak 
enough that we cannot exclude the possibility that 
the TeV, if not the GeV, emission is due to inverse- 
Compton scattering [ Aliu el aL l l2013t . On the other 
hand, it is also plausible that the VHE 7-ray emis- 
sion is produced by interaction of accelerated nuclei 
with the H i shell surrounding the remnant. Es- 
timates of the target material density required for 
accelerated nuclei to produce the obs erved VHE 7- 
ray flux, based on iDrurv et al. Il994 [ give a range 



of densities 1.0 — 5.5 cm" -3 |Aliu el al.l 120131 consis- 
tent ^ Gosachinskiil 2001 1 's estimates of the gas density 
within the Hi shell. However, it must be noted that 
the shock velocities inferred from the optical and X- 
ray data are too low for the forward shock to be cur- 
rently accelerating particles to TeV energies. If the 
VHE 7-ray emission is hadronic, it is likely due to 
particles accelerated when the remnant was younger 
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Figure 4: Top Panel — ASCA X-ray spectrum of the 
region of enhanced X-ray emission coincident with 
VER J2019+407 shown in Figure [3 The stepped line 
shows the fit of a Raymond-Smith thermal plasma model 
with parameters as given in the text. Bottom Panel — 
residuals from the best-fit model. 



that are now interacting with the shell. 

The relationship of the cocoon of freshly-accelerated 
cosmic rays detected by Fermi to SNR G78. 2-1-2.1 and 
VER J2019-I-407 also remains unclear. It is possible 
SNR G78. 2+2.1 either has injected or is injecting ac- 
celerated particles into the cocoon. However, while 
they are shown for reference , we caution against us 
ing the cocoon contours from I Ackermann et al.l 2011 



to judge the relationship of the cocoon to the VHE 7- 
ray emission, since they are derived from an analysis 
where IFGL J2020. 0-1-4049, which is no longer consid- 
ered an independent source, was included as part of 
the background model. It should also be noted that 
the VHE 7-ray excess map in this paper was made 
with the ring-background estimation method, which is 
ill-suited to detecting a large-scale 4 square degree) 
region of 7-ray emission such as the cocoon. Therefore 
the VHE 7-ray maps shown here cannot be used to set 
a meaningful upper limit on cocoon emission above 
300 GeV. A conclusion determination of the relation- 



ship between SNR G78. 2+2.1 and VER J2019+407 
will have to await further data, analyzed with more 
sophisticated analysis techniques. 
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